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COMPARATIVE  INVESTIGATION  OF  HIGH-LATITUDE  IONOSPHERIC 
STRUCTURE  AND  EFFECTS  NEAR  SOLAR  MAXIMUM 


L  INTRODUCTION  AND  OBJECTIVES 

Northwest  Research  Associates  (NWRA)  has  recorded  a  database  containing  measure¬ 
ments  of  the  intensities  of  two  L-band  signals  transmitted  by  satellites  of  the  Global  Positioning 
System  (GPS)  and  the  differential  carrier  phase  (DCP)  and  differential  group  delay  (DGD) 
between  them.  Data  collection  was  carried  out  during  the  period  December  1990  -  December 
1992,  which  immediately  followed  the  recent  solar  maximum  This  database  contains  information 
collected  at  Thule  AB,  Greenland;  Lerwick  Observatory,  Shetland  Island,  UK;  and  Hanscom 
AFB,  MA,  and  represents  a  range  of  latitude  and  longitude  regions.  Software  has  been  developed 
that  processes  these  data  and  produces  measurements  of  trans-ionospheric  absolute  total  electron 
content  (TEC)  and  L-band  scintillation.  This  database  contains  morphology  of  TEC  in  the  polar 
cap  and  in  the  region  of  the  BMEWS  SITE-ID  radar  at  Fylingdales,  UK,  near  solar  maximum. 

The  purpose  of  this  effort  was  to  complete  a  database  sufficient  to  characterize  TEC 
variations  and  scintillation  at  high  latitudes  near  solar  maximum  and  to  document  the  behavioral 
trends  therein.  One  specific  objective  of  the  project  was  to  complete  a  database  that  allows 
improved  characterization  of  TEC  variations  and  scintillation  at  high  latitudes  near  solar  maxi¬ 
mum.  These  data  are  suitable  for  analyses  aimed  at  both  Air  Force  systems  at  transition  and 
higher  latitudes  and  understanding  of  ionospheric  behaviors  that  affect  those  operations.  The 
second  objective  was  to  advance  the  behavioral  description  and  understanding  of  plasma-density 
structures  within  geomagnetic  latitude  regimes  identified  with  the  main  ionospheric  trough,  the 
auroral  oval,  and  the  polar  cap  during  years  near  the  peak  of  the  eleven-year  cycle  in  solar  activ¬ 
ity.  This  will  form  a  basis  for  operational  and,  potentially,  design  decisions  regarding  surveillance 
radars,  with  possible  applications  to  transionospheric  communications  systems. 

To  accomplish  these  goals,  NWRA  has  maintained  and  operated  (or  directed  the  operation 
of)  receiving  equipment  at  Thule  AB,  Greenland;  Lerwick  Observatory,  Shetland  Island,  UK;  and 
Hanscom  AFB,  MA.  The  received  TEC  and  scintillation  data  have  been  stored  on  various  mag¬ 
netic  media.  Software  was  developed  to  reduce  and  analyze  the  data  to  advance  the  under¬ 
standing  of  the  nature  and  development  of  ionospheric  structures  within  the  regions  monitored. 
Software  was  also  developed  to  apply  the  database  to  characterizing  the  transionospheric  radio 
channel. 


a  DATA  COLLECTION  AND  PROCESSING 

During  the  22-month  study  period,  NWRA  maintained  and  operated  (or  directed  the 
operation  of)  data-collection  equipment  at  Thule  AB,  Greenland;  Lerwick  Observatory,  Shetland 
Island,  UK;  and  Hanscom  AFB,  MA. 
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A.  Description  of  Instrumentation 


The  government  furnished  equipment  (GFE)  operated  and  serviced  by  NWRA  consists  of 
Stanford  Telecommunications,  Inc.  (STEL)  STEL-5010  single-channel  GPS  receivers,  a  Texas 
instrument  (TT)  TI-4100  four-channel  GPS  receiver,  Magnavox  (MX)  MX-1502  Transit  re¬ 
ceivers,  and  a  National  Institute  of  Standards  and  Technology  Ionospheric  Measurement  System 
(NIMS)  code-free,  multi-channel  GPS  receiver  manufactured  by  Atmospheric  Instrumentation 
Research,  Inc.  (AIR). 

The  STEL-5010  located  at  Shetland  Island  experienced  many  failures  during  this  study 
period.  Many  resources,  both  financial  and  manpower,  were  expended  to  maintain  the  site  in  an 
operating  condition.  This  impacted  the  quantity  and  quality  of  data  that  were  collected  at  the  site, 
as  well  as  the  amount  of  data  from  all  sites  that  were  processed.  Emphasis  was  placed  on  die 
reduction  of  TI-4100  four-channel  GPS  Shetland  Island  data,  which  are  applicable  to  the 
BMEWS  radar  at  Fylingdales,  UK. 

The  STEL-5010  single  channel  GPS  receiver  tracks  a  single  satellite  for  up  to  six  hours. 
It  measures  the  DGD  between  LI  (1575.4  Mhz)  and  L2  (1227.6  Mhz),  the  DCP  between  those 
two  signals,  and  the  received  intensities  of  the  LI  and  L2  signals  from  the  older.  Block- 1  GPS 
satellites  only.  It  cannot  receive  the  Block-2  satellites,  whose  L2  signals  are  encoded  so  the  GPS 
selective  availability  function  can  be  invoked.  The  data  were  recorded  at  twenty  samples  per 
second  as  binary  files. 

The  TI-4100  four-channel  receiver  collected  data  simultaneously  from  four  different  GPS 
satellites,  both  Block  1  and  Block  2.  DGD  and  DCP  were  recorded  at  one  sample  per  second. 
However,  the  DCP  data  experienced  many  drop  outs,  making  these  data  unsuitable  to  phase  aver¬ 
aging  as  a  method  of  referencing.  The  DGD  data  were  essentially  continuous  and,  with  proper 
antenna  siting,  contained  minimal  multipath.  These  DGD  data  were  used  to  generate  equivalent 
vertical  TEC.  Intensities  of  the  LI  and  L2  signal  are  not  directly  measured,  therefore  scintillation 
indices,  S4,  cannot  be  derived  from  the  TI-4100  data. 

The  MX-1502  Transit  receivers,  located  at  Shetland  Island,  UK,  and  Hanscom  AFB,  MA, 
provide  latitudinal  relative  TEC  from  north-south  traveling  satellites  of  the  Navy  Navigation 
Satellite  System  (NNSS).  The  transmitted  signals  from  these  satellites  are  at  150  Mhz  and  400 
Mhz.  The  sampling  rate  of  this  receiver  is  not  fast  enough  to  yield  scintillation  information.  Each 
satellite  pass  lasts  approximately  fifteen  minutes.  In  order  to  determine  absolute  TEC  levels  from 
the  Transit  measurements,  which  only  provide  relative  TEC  along  the  satellite  track,  comple¬ 
mentary  absolute  TEC  data  are  required  to  calibrate  the  relative  TEC  measurements.  Absolute 
TEC  information  from  the  co-located  GPS  receivers  provides  such  data. 

The  NIMS  receiver  is  a  code-free  GPS  receiver  whose  output  is  fifteen-minute  averaged 
values  of  relative  TEC.  This  receiver  was  acquired  in  April  1992  but  required  repair  and  was 
returned  to  AIR,  Inc.  Mr.  Charley  Andreasen  (NWRA)  identified  several  defects,  all  of  which 
were  confirmed  and  corrected  by  the  manufacturer.  The  unit  was  installed  at  Hanscom  AFB  upon 
its  return  in  May  1992.  In  a  letter  dated  27  July  1992,  Mr.  David  B.  Call,  president  of  AIR,  Inc., 
reported  that  this  receiver  fails  to  meet  the  specified  accuracy  levels.  The  problem  manifests  itself 
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in  rapid  fluctuations  in  measured  TEC  that  are  due  to  the  antenna  and  choke-ring  assembly,  not 
the  ionosphere.  AIR  intends  to  correct  this  problem  with  a  redesigned  antenna.  However,  the 
data  from  this  unit  are  necessarily  of  doubtful  quality  due  to  this  design  problem. 

B.  Instrumentation  and  Operation  at  Each  Site 

THULE  AB 

One  STEL-5010  single-channel  GPS  receiver  was  operated  at  the  Thule  AB  site.  The 
receiver  reliably  collected  data  throughout  most  of  the  study  period.  Based  on  sample  data  that 
have  been  processed,  the  quality  appears  to  be  good.  Appendix  4,  Near  Solar  Maximum  Data¬ 
base,  may  be  referenced  to  determine  coverage  for  a  specific  date.  Information  pertaining  to  the 
maintenance,  equipment  calibration,  problems,  and  repairs  at  this  site  are  detailed  in  Appendix  1. 

SHETLAND  ISLAND 

At  the  Shetland  Island  site,  NWRA  directed  the  operation  of  three  receivers:  an  STEL- 
5010  single-channel  GPS  receiver,  a  TI-4100  four-channel  GPS  receiver,  and  an  MX- 1502 
Transit  receiver.  In  addition,  a  Transit  receiver  operated  by  the  University  College  of  Wales 
(UCW),  Aberystwyth,  UK,  collected  data  at  this  site.  The  STEL-5010  experienced  many  failures 
during  the  study  period,  and  reliable  TEC  and  scintillation  data  are  available  only  from  the 
November  1991 -June  1992  period.  The  TI-4100  operated  reliably  from  December  1990  -  Feb¬ 
ruary  1992,  when  it  was  sent,  at  Phillips  Laboratory's  request,  to  Applied  Research  Laboratories, 
University  of  Texas  at  Austin  (ARL-UT),  to  be  calibrated  and  reprogrammed.  Absolute  TEC 
data  quality  for  this  period  is  good.  Data  collected  by  the  MX- 1502  Transit  receiver  are  available 
from  most  of  the  November  1991 -June  1992  period,  and  are  of  good  quality.  UCW  Transit  data 
arc  available  from  most  of  the  study  period.  Appendix  4  may  be  referenced  to  determine  cover¬ 
age  for  a  specific  date.  Information  pertaining  to  the  maintenance,  equipment  calibration,  prob¬ 
lems,  and  repairs  at  this  site  are  detailed  in  Appendix  1. 

UCW  references  their  Shetland  Island  Transit  data  to  the  GPS  data,  both  single-  and  four- 
channel,  collected  at  Shetland  Island  by  NWRA.  The  programs  that  plot  the  GPS  data  generate 
tables  of  time,  TEC,  latitude,  and  longitude,  at  five-minute  intervals,  which  are  used  by  UCW  as 
reference. 

HANSCOM  AFB 

Three  receivers  were  operated  at  the  Han  scorn  AFB  site:  an  STEL-5010  single-channel 
receiver,  an  MX  1502  Transit  receiver,  and  an  NIMS  code-free,  multichannel  GPS  receiver.  The 
STEL-5010  collected  data  during  the  period  December  1 990-October  1991.  Samples  of  data  that 
have  been  processed  indicate  a  problem  with  the  DGD  channel  dropping  out  intermittently, 
making  determination  of  equivalent  vertical  TEC  from  phase-averaged  DCP  impossible  during 
those  dropouts.  The  Transit  receiver  was  installed  in  July  1991  and  collected  data  through  some 
of  the  remaining  study  period.  The  NIMS  receiver  was  in  operation  from  May  1992-December 
1992,  but,  as  discussed  previously,  antenna  and  choke-ring  assembly  design  flaws  resulted  in  poor 
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data  quality.  Information  pertaining  to  the  maintenance,  equipment  calibration,  problems,  and 
repairs  at  this  site  are  detailed  in  Appendix  1. 

In  the  Han  scorn  longitude  sector,  absolute  TEC  data  for  referencing  the  Transit  measure¬ 
ments  are  available  from  sources  other  than  the  co- located  receivers.  From  January  1991  through 
June  1992,  data  from  a  geostationary  satellite  were  collected  using  a  polarimeter  located  at 
Hamilton,  MA.  These  data  can  be  requested  through  Mr.  Jack  Klobuchar  at  AFPL/GPIM, 
Han  scorn  AFB,  MA.  Ionosonde  data  collected  at  Millstone  Hill,  MA,  are  available  in  several  dif¬ 
ferent  formats  from  World  Data  Center,  Boulder,  CO,  and  could  also  be  used  as  reference. 
(Ionosonde  data,  collected  during  the  contract  period,  with  approximately  90  %  yearly  coverage, 
are  also  available  from  Qanac,  Greenland,  and  Narssarssuaq,  Greenland.) 

C.  Absolute  Total  Electro  *  Content  Data 

Originally,  it  was  planned  to  generate  mean  absolute  TEC  as  one  of  the  summary  parame¬ 
ters  contained  in  the  processed  database.  Upon  further  consideration,  it  was  decided  that  that 
value  is  not  directly  useful,  due  to  the  GPS  satellites'  geographic  and  local  time  spread.  It  was 
decided  instead  to  analyze  the  data  by  separating  them  by  latitude  regions  and  plotting  diem.  In 
this  way,  the  behavior  associated  with  the  different  latitude  regions  (trough  edge/auroral,  or 
auroral/polar  cap)  can  readily  be  seen.  In  addition,  the  latitude-separated  data  may  be  plotted 
according  to  Ionospheric  Penetration  Point  (IPP)  adjusted  time.  To  determine  the  IPP-adjusted 
time,  the  latitude  and  longitude  of  the  IPP  are  first  calculated.  A  time  offset  is  added  to  the 
Universal  Tune  Coordinated  (UTC)  time  to  reflect  the  time  shift  associated  with  the  location  of 
the  IPP.  Using  this  option,  the  TEC  curves  from  two  different  satellites  whose  paths  traverse  the 
same  region  of  the  ionosphere  will  tend  to  overlap.  This  results  in  much  more  readily  analyzed 
24-hour  TEC  profiles.  Figure  1  shows  four  channels’  (trackers)  data  from  the  TI-4100  plotted 
individually  with  latitude  and  longitude  of  the  IPP.  A  measure  of  variability  of  TEC  with  look 
direction  is  evident  For  comparison.  Figure  2  shows  the  same  data  that  have  been  separated  by 
latitude  region  and  plotted  versus  IPP-adjusted  time.  Another  comparison  of  UTC-plotted  data 
and  IPP-adjusted  time  plotted  data  is  given  in  Figure  3.  STEL-5010  single-channel  data  from 
Days  16  through  23, 1992  are  overplotted  in  both  UTC-  and  EPP-adjusted  time. 

Measurements  of  DGD,  from  which  absolute  TEC  can  be  determined,  are  very  susceptible 
to  contamination  by  multipath.  The  DCP  yields  a  relative  measure  of  TEC,  containing  a  2an 
ambiguity.  By  referencing  the  DCP  to  the  DGD  through  the  phase  averaging  process,  equivalent 
vertical  TEC  can  be  determined  from  DCP  data. 

TEC  morphology  for  the  TI-4100  was  derived  from  DGD  alone,  due  to  the  intensive  labor 
input  that  would  have  been  required  to  perform  phase  linking  of  the  discontinuous  DCP  data  to 
support  the  phase-averaging  process.  Multipath  is  present  in  this  DGD  data,  but  due  to  the 
antenna  having  been  located  in  an  open  field,  its  magnitude  is  about  ±  2  to  3  TEC  units.  The 
STEL-5010  data,  on  the  other  hand,  had  continuous  DCP  data  on  most  of  the  passes,  and  this 
was  used  to  eliminate  multipath  from  the  absolute  TEC  data. 

The  program  that  plots  the  data  according  to  latitude  could  be  further  modified  to  gen¬ 
erate  a  plot  of  the  mean  diurnal  TEC  curve  over  some  desired  interval,  for  example,  one  week. 
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Figure  2.  TI-4100  Data  separated  by  latitude. 
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Figure  3.  STEL-5010  Data  plotted  vs  UTC  and  IPP  adjusted  time. 
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To  do  this,  the  TEC  values  at  each  time  for  each  data  day  could  be  averaged  together  in  an  array 
during  overplotting,  and  that  array  could  then  be  plotted  as  a  mean  diurnal  TEC  curve  after  the 
requested  several-day  overplot  was  produced. 

D.  Approach  for  Statistical  Assessment  of  TEC  Variation  and  Scintillation 

During  periods  of  high  planetary  Kp  index,  the  occurrence  of  scintillation  and  the  mag¬ 
nitude  of  TEC  variations  are  increased,  making  these  periods  especially  valuable  for  study.  First- 
order  .mstical  measures,  such  as  standard  deviation  of  relative  TEC  variation  and  L-band 
intensity  scintillation,  are  quantities  of  interest  Single-channel  data  from  the  STEL-5010  systems 
are  best  suited  to  an  analysis  of  these  quantities  since  it  alone  provides  scintillation  data  and  since 
the  four-channel  data,  as  noted  before,  suffers  breaks  in  phase  continuity  which  would  be  very 
labor  intensive  to  correct  Equivalent  vertical  TEC  that  has  been  derived  from  phase-averaged 
and  slant-factor  corrected  DCP  could  be  detrended,  and  the  standard  deviation  could  then  be 
calculated.  It  is  necessary  to  use  phase-averaged,  slant-factor  corrected  DCP  data  since  noise 
also  scales  up  at  low  elevation. 

E.  Software  Tools  Developed  for  Data  Analysis 

Appendix  3  contains  a  description  of  the  software  that  was  developed  to  process  these 
data  into  a  form  suitable  for  analyses  by  both  Air  Force  systems  at  transition  and  higher  latitudes, 
as  well  as  improved  understanding  of  ionospheric  behaviors  that  affect  those  operations.  These 
programs  can  be  used  with  any  GPS  data,  in  the  same  manner  as  they  are  used  with  the  NWRA 
database:  to  test,  quality  check,  reduce,  analyze,  and  display  the  data. 

The  software  developed  during  this  effort  reduces  the  raw  data  to  files  from  which  24- 
hour  plots,  several-day  overplots,  pass-file  plots,  latitude-separated  plots,  and  IPP-adjusted  time 
plots  of  equivalent  vertical  TEC  can  be  produced.  S4  at  LI  and  L2  can  be  plotted  as  24-hour 
plots  or  individual  pass-file  plots.  A  summary  of  the  STEL-5010  data  in  tabular  format  can  also 
be  generated.  This  summary,  shown  in  Figure  4,  contains  universal  time  of  observation  in 
seconds,  amplitude  scintillation  index  (S4)  at  LI  and  L2,  slant  path  TEC,  RMS  dispersive  phase 
from  DCP  measurements,  azimuth  and  elevation  angle  of  the  satellite  from  the  observing  station, 
and  geomagnetic  latitude,  longitude,  time,  and  penetration  angle  at  the  3 50- km  IPP.  These 
programs  can  be  used  with  any  GPS  data,  such  as  those  collected  by  the  planned  Trans- 
Ionosphcric  Sensing  System  (TISS),  in  the  same  manner  as  they  are  used  with  the  NWRA  data¬ 
base:  to  test,  quality  check,  reduce,  analyze,  and  display  the  data. 

A  complete  description  of  the  available  software  to  process  the  near  solar  maximum  data¬ 
base  resulting  from  this  effort  is  given  in  Appendix  3.  These  programs  can  be  used  with  any  GPS 
data,  such  as  that  collected  by  the  planned  TISS. 

F.  Database 

The  near-solar-maximum  database  in  Appendix  4  details  the  availability  of  data  according 
to  date,  site,  and  receiver  type,  and  indicates  the  level  to  which  the  data  have  been  processed  and 
the  quality  of  that  processed  data.  Each  of  the  25  pages  shows  one  month  broken  down  into 
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Figure  4.  Sample  of  summary  parameters  -  output  of  GPSPCNW3.FOR. 


days.  Specified  by  site,  the  listing  includes  a  column  for  each  receiver  type.  An  X  indicates  that 
raw  data  from  that  site,  on  that  day,  collected  by  that  receiver,  is  available.  An  O  appears  if  no 
data  exist  In  a  separate  column,  Y  for  yes  and  N  for  no  indicates  if  the  data  have  been  processed. 
The  last  column  shows  the  quality  of  the  data  that  have  been  processed:  0  if  data  not  processed;  1 
for  good;  2  for  questionable  but  possibly  recoverable;  and  3  for  bad. 


EL  RESULTS 

A.  TEC  Results:  Shetland  Island 

The  data  collected  at  Shedand  Island  by  the  TI-4100  four  channel  GPS  receiver  yielded  a 
great  deal  of  data  for  studying  TEC  only.  The  TEC  morphology  of  Shetland  Island  from  March 
1991  to  February  1992  was  determined  by  overplotting  approximately  one-week  sets  of  equiva¬ 
lent  vertical  TEC  derived  from  DGD.  These  are  presented  in  Figures  5-16.  As  noted  before, 
use  of  DGD  alone  was  necessary  due  to  the  intensive  labor  effort  that  would  have  been  required 
to  perform  phase  linking  on  the  DCP  data  in  order  to  support  the  phase-averaging  process.  How¬ 
ever,  since  the  antenna  was  located  in  an  open  field,  these  curves  contain  multipath  noise  of  about 
±  2  to  3  TEC  units.  For  cases  of  special  interest  (a  few  days),  it  is  possible  and  reasonable  to  use 
the  PLTREC7C.FOR  program,  link  the  DCP,  and  perform  phase  averaging.  This  would  effec¬ 
tively  eliminate  multipath  noise  and  the  data  would  look  like  the  TEC  from  STEL-5010  single 
channel  GPS  data. 

A  review  of  the  year  of  weekly  Shetland  Island  TI-4100  four-channel  overplots  in  Figures 
5  -16  shows  that  the  TEC  profiles  exhibit  variations  dependent  on  time  of  day,  season  of  the  year, 
latitude,  and  Kp.  The  24-hour  plots  of  TEC  show  the  typical  cosine-like  diurnal  behavior  of  the 
mid-latitude  ionosphere  starting  around  0900  hours  and  lasting  until  1800  with  the  peak  occurring 
at  approximately  1200. 

The  data  plotted  begin  in  March  1991,  where  TEC  levels  were  usually  quite  high  during 
the  day  and  there  was  significant  variability,  which  is  consistent  with  the  large  magnetic  storm  that 
occurred  that  month.  April  transitions  from  the  high  daytime  TEC  levels  of  March  to  the  nearly 
"flat"  diurnal  behavior  that  tends  to  be  seen  during  the  summer  months  in  the  mid-latitude  iono¬ 
sphere.  This  flat  behavior  was  seen  to  persist  until  September,  which  transitions  to  a  clear  post¬ 
midday  peak  and  higher  levels,  though  lower  than  was  seen  in  March.  The  month  of  October 
tends  to  exhibit  high  TEC  activity  levels,  with  the  last  week  of  the  month  showing  high  variability 
due  to  a  magnetic  storm.  November  shows  much  disturbance,  having  high,  medium,  and  low 
diurnal  peaks.  Kp  levels  were  high  during  much  of  this  month,  with  a  large  storm  on  Day  312. 
December  shows  very  stable  classic  diurnal  TEC  variation  throughout  the  month,  with  moderate 
maximum  levels.  This  again  is  consistent  with  the  lower  Kp  levels  seen  in  December.  January 
and  February  1992  show  more  variability  and  higher  TEC  levels.  Very  high  TEC  levels  were  seen 
associated  with  the  storm  period  of  02  -  03  and  08  February. 

During  periods  of  increased  solar  activity  corresponding  to  high  Kp,  rapid  variations  in 
TEC  levels  and  higher  scintillation  values  are  seen.  Also,  it  is  not  uncommon  to  see  the  signature 
of  the  trough  in  the  24-hour  TEC  profiles.  This  is  particularly  evident  in  the  01  -  08  February 
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Figure  5.  Ladtudinally  separated  TT-4100  data  -  March  1991. 


I 


o 

in 

u» 

i- 

► 

O 

u 

III 

I- 

I- 

n-: 

tij 

> 

> 

CJ 

in 


°0  3  6  9  12  IS  18  21  24 

LOCAL  TIME  AT  IPP 


-J - 1 - 1 - 1 - 1 - 1 - 1- 


I.JAY  16791 


SHETLAND  - 
FOUR  TRACKS 
SOUTH 


— i — i — i — i — i — i — i — 


DRY  17B91 


SI  C  l  LAND  - 
FOUR  TRACKS 
SOUTH 


0  3  G  9  12  IS  in  21  :m 

LOCAL  TIME  AT  IPP 


Figure  8.  Ladtudinally  separated  TI-4100  data  -  June  1991. 
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Figure  11.  Latitudinally  separated  TI-4100  data  *  September  1991. 
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Figure  12.  Lalitudinally  separated  TI-4100  data  -  October  1991. 
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Figure  13.  Latitudinally  separated  TI-4100  data  -  November  1991. 
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Figure  14.  Latitudinally  separated  TI-4100data  •  December  1991. 


Figure  15.  Latitudinally  separated  TI-4100data  -  January  1992. 


22 


overplot  shown  in  Figure  16  which  includes  the  02  -  03  and  08  February  storms.  The  trough,  a 
region  of  lower  TEC  density  in  the  ionosphere,  is  typically  located  poleward  from  the 
midlatitudes.  It  is  seen  to  move  equatorward  during  periods  of  high  Kp,  and  its  signature  is 
characterized  by  steepened  gradients  of  TEC. 

B.  TEC  Results:  Other  Sites 

Similar  overplots  from  single  channel  Shetland  Island  data  can  be,  and  in  some  cases  have 
been,  generated.  These  plots  show  the  same  behaviors  as  were  seen  in  the  Tl-4100  data 

TEC  morphology  from  Thule  AB  single-channel  data,  displaying  behavior  of  the  iono¬ 
sphere  in  the  polar  cap  region,  can  be  produced  from  the  database  using  the  software  tools  devel¬ 
oped  under  this  effort 

Samples  of  the  single-channel  data  from  the  Han  scorn  site  indicate  that  there  was  a  prob¬ 
lem  with  the  DGD  data  channel  The  extent  of  this  problem  will  not  be  known  until  more  data  are 
analyzed.  However,  should  the  data  quality  for  some  of  the  measurements  be  found  to  be  good, 
similar  overplots  from  Han  scorn  can  be  produced  to  provide  insight  into  ionospheric  behavior  in 
that  region. 

C.  TEC  Results:  Storm  Time  Examples 

Data  collected  by  the  TI-4100  four  channel  GPS  receiver  from  two  magnetically  active 
periods  in  November  1991  are  shown  in  Figures  17  and  18.  A  chart  of  Kp  level  corresponding  to 
these  days  is  found  in  Figure  19.  Days  301-307,1991  (28  October  -  3  November,  1991)  and  days 
311,  312,  and  313,  1991  (7-9  November  1991)  have  been  plotted  according  to  north,  overhead, 
and  south  latitude  regions  relative  to  the  observing  station's  location.  Plots  of  the  magnetically 
quiet  day  311  show  the  typical  cosine-like  midlatitude  behavior.  At  0100  to  0500  hours,  the 
trough-wall  signature  is  seen.  There  is  a  sudden  commencement  on  day  312,  and  Kp  levels  rise  to 
a  high  of  nine  over  the  next  24  hours.  The  trough-wall  signature,  characterized  by  steep  gradients 
in  TEC,  is  evident  in  these  plots.  Day  313  shows  a  relatively  flat  TEC  profile  in  the  latter  half  of 
the  day  as  the  magnetic  storm  ends  and  Kp  levels  return  to  zero.  The  technique  of  plotting 
latitudinally  separated  data  allows  viewing  of  the  equatorward  movement  of  the  trough.  In  Figure 
17,  the  trough-wall  signature  is  seen  at  approximately  1500  hours  in  the  overhead  direction,  and 
appears  an  hour  later  in  the  southerly  data.  Increased  levels  of  TEC  and  the  appearance  of  the 
trough  wall  are  also  seen  in  the  overhead  plots  in  Figure  18.  These  plots  show  distinct  features  in 
each  latitude  region.  The  overhead  plot  of  days  306  and  307,  as  the  storm  is  subsiding,  shows  die 
return  to  a  more  normal  cosine-like  shape  and  lower  TEC  values. 

D.  Scintillation  Data 

Single-channel  data  are  available  from  all  three  sites  concurrently,  for  a  number  of  days  in 
1991  during  periods  of  high  planetary  Kp.  These  times  are  March  24,  25,  1991;  June  5,  6,  9,  10, 
17,  1991;  and  October  28,  29,  31,  1991.  A  comparative  analysis  of  the  TEC  profiles  and  scin¬ 
tillation  between  the  two  longitude  sectors  could  be  done  using  these  data. 
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Figure  17.  TT-4100  data  -November  7-9, 1991  storm  period.  LOCAL  TIME  OT  ipp 
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Figure  18.  TI-4100  data  -  October  28  -  November  4, 1991 


An  analysis  was  performed  on  measurements  taken  with  the  S TEL- 50 10  receiver  at 
Lerwick  Observatory,  Shetland  Island,  UK,  on  Day  327,  1991  to  Day  60,  1992.  Raw  data  were 
reduced  using  the  GPSPCNW3POR  program  to  produce  propagation  parameters  associated  with 
the  transionospheric  propagation  path  and  consisting  of  amplitude  scintillation  index  (S4)  at  LI 
and  L2,  TEC  over  the  slant  path  (in  units  of  1016  e/m2),  and  the  RMS  variation  in  phase  (also  in 
equivalent  TEC  units).  Geometry  information  stored  with  these  parameters  to  facilitate  data 
analysis  consists  of  the  elevation  angle  and  azimuth  of  the  satellite  as  seen  from  the  ground 
station,  the  geomagnetic  latitude  and  longitude  coordinates  of  the  350-km  altitude  IPP,  the 
incidence  angle  at  the  IPP,  and  geomagnetic  time  at  the  IPP.  Geomagnetic  coordinates  used  refer 
to  the  magnetic  apex  (Secan,  1987;  VanZandt  et  al,  1972). 

A  limited  study  was  performed  that  concentrated  on  the  review  and  interpretation  of  avail¬ 
able  measurements  during  active  periods.  Initially  the  database  was  scanned  and  periods  were 
identified  when  S4(L2)  was  greater  than  0.20  for  elevation  angles  of  30  deg  or  greater.  A  sum¬ 
mary  of  episodes  of  enhanced  scintillation  and  various  measurement  characteristics  is  listed  in 
Table  1.  As  shown,  amplitude  scintillation  index  at  L2  ranged  from  0.261  to  a  maximum  0.558 
and,  in  all  but  one  instance,  corresponded  to  high  levels  of  planetary  Kp  index,  which  ranged  from 
2+  to  8.  The  geomagnetic  latitude  at  the  IPP  intersection  was  mostly  in  the  sub-auroral  zone  ’uid 
varied  from  55  to  62  deg.  The  measurements  at  Lerwick  Observatory,  Shetland  Island,  UK, 
normally  show  quiet  scintillation  conditions  consistent  with  observations  south  of  the  auroral-zone 
scintillation  boundary.  During  periods  of  activity,  in  most  instances  coincident  with  times  of  high 
planetary  Kp  index,  the  observations  are  consistent  with  the  equatorward  movement  of  the 
irregularity  boundary. 


TIME  GEOMAGNETIC  MAX.%  TIME  GEOMAGNETIC 


1992 

DAY 

SAT 

MAX  (UT) 
S4(L2)  (hours) 

TIME 

(hours) 

LAT. 

(deg) 

CHANGE 

TEC(V) 

(UT) 

(hours) 

TIME 

(hours) 

LAT. 

(deg) 

Kp 

16 

3 

0.405 

2320 

23.82 

57.89 

39.7 

23.08 

23.68 

57.67 

5+ 

33 

11 

0.292 

2.78 

3.35 

59.41 

6.8 

2.63 

3.19 

59.49 

6 

33 

3 

0.307 

23.57 

0.49 

59.55 

7.1 

23.52 

0.43 

59.53 

6 

34 

11 

0.314 

3.68 

4.32 

58.48 

13.6 

3.68 

4.31 

58.48 

7 

57 

3 

0.558 

19.80 

20.43 

56.59 

11.6 

19.78 

20.41 

56.54 

8 

60 

13 

0.433 

4.93 

5.38 

63.57 

3.0 

5.45 

6.01 

61.66 

2+ 

60 

3 

0261 

20.10 

20.78 

57.67 

6.4 

20.11 

20.78 

57.68 

5+ 

Table  1 .  Summary  of  Measurements  During  Active  Periods  at  Shetland. 
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In  each  of  the  seven  instances  where  enhanced  scintillation  was  observed,  the  measure¬ 
ments  are  summarized  using  six  graphic  presentations,  which  are  presented  in  Figures  20  -  26. 
The  upper  left-hand  figure  in  each  presentation  displays  the  azimuth/elevation  track  of  the  satellite 
during  the  period  of  activity.  The  upper  right-hand  figure  shows  the  geomagnetic  latitude  as  a 
function  of  geomagnetic  time.  The  middle  plots  show  the  equivalent  vertical  TEC  (left  side)  and 
apparent  fractional  variation  in  slant  TEC  (right  side)  as  a  function  of  universal  time.  The  bottom 
figures  show  plots  of  S4  index  at  LI  and  L2  (left  side)  and  geomagnetic  latitude  (right  side)  as  a 
function  of  universal  time. 

With  one  exception,  the  scintillation  analyzed  occurred  between  four  hours  before  and 
after  midnight  The  isolated  activity  event  occured  at  five  hours  after  midnight,  and  it  differed 
also  in  that  it  was  accompanied  by  only  a  2+  level  of  Kp.  In  all  the  events,  the  S4  level  was  higher 
at  the  lower  frequency  (L2),  as  expected  from  propagation  theory. 

The  reported  fractional  change  in  slant  TEC  is  available  from  the  preliminary  processing  at 
a  20  seconds  per  measurement  interval  and  is  the  RMS  measured  TEC  obtained  at  the  20  samples 
per  second  (SPS)  rate  over  a  one-minute  measurement  interval  normalized  by  the  average  TEC 
As  defined  here,  the  measured  apparent  fractional  variation  in  TEC  can  be  either  a  true  variation 
in  TEC  (a  refraction  process)  or  indicative  of  a  diffraction  process  or  it  can  be  a  combination  of 
both  refaction  and  diffraction. 

For  example,  in  reviewing  the  graphical  data  covering  the  enhanced  activity  period  on  Day 
16,  the  sharp  change  in  vertical  TEC  near  2300  hours  appears  to  be  clearly  related  to  large-scale 
refraction  ionospheric  effects.  On  the  other  hand,  the  rapid  variation  in  the  apparent  factional 
TEC  shown  in  the  next  group  of  plots  covering  Day  33  shows  rapid  smaller  fractional  changes  in 
TEC  that  appear  to  be  more  characteristic  of  a  diffraction  process. 

E.  Examples  of  TEC  Data  vs  Model  Prediction 

Since  the  character  of  the  ionosphere  in  the  Shetland  Island  region  is  seen  to  vary,  it  is  of 
interest  to  consider  the  accuracy  of  the  model  that  is  used  to  correct  the  range  error  of  radar 
measurements.  A  statistical  comparison  between  the  USAF  Environmental  Tactical  Application 
Center  (ETAC)  ionospheric  model,  based  on  the  Bent  ionospheric  model  (Llewellyn  and  Bent 
1973)  and  actual  measured  values  of  TEC  is  accomplished  using  the  Transit  dialysis  software 
package  developed  by  RDP,  Inc.,  Waltham,  MA.  The  statistical  analysis  is  displayed  as  scatter 
plots  of  Bent-model  predicted  values  of  TEC  versus  measured  values  of  TEC,  and  as  histograms 
depicting  the  relative  percent  error  between  those  two  quantities.  The  data  can  be  analyzed 
according  to  a  specific  elevation  range,  azimuth  range,  or  time  period. 

The  Transit  analysis  software  was  originally  developed  to  work  with  Transit  satellite  data 
collected  at  Lerwick  Observatory,  Shetland  Island,  UK,  by  the  UCW  and  referenced  to  the  GPS 
data  collected  at  Lerwick  Observatory,  Shetland  Island,  UK,  by  NWRA.  It  has  been  modified  to 
accept  also  GPS  data  from  Lerwick  Observatory,  Shetland  Island,  UK,  directly,  as  in  the  exam¬ 
ples  of  scatter  plots  and  histograms  shown  in  Figures  27  -  32.  One  week  of  Shetland  GPS  data 
from  the  TI-4100  four  channel  receiver  from  December  1991  are  analyzed. 
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Figure  21.  Graphical  analysis  of  S TEL- 5010  Shetland  data  -  Day  33. . 
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Figure  25.  Graphical  analysis  of  STE^-5010  Shetland  data  -  Day  60. 
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Figure  26.  Graphical  analysis  of  STEL-5010  Shetland  data  -  Day  60. 
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Figure  29.  Comparison 
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Figure  31.  Comparison  of  Bent  model  predictions 
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Figure  32.  Comparison  of  Bom  model  predictions  to  measured  TEC  -  Shetland  December  199ITI-4100dau-  histogram 


F.  Data  Use  in  Scintillation  Model 

A  potential  use  for  measurements  of  the  effects  of  the  ionosphere  on  GPS  signals  is  to  use 
observations  of  ionospheric  scintillation,  both  intensity  and  phase,  to  update  a  global  model  of 
ionospheric  scintillation  in  near  real  time.  The  target  model  is  the  WBMOD  ionospheric  scintil¬ 
lation  model  presently  in  operational  use  by  the  USAF  Air  Weather  Service  (AWS).  In  order  to 
use  GPS-derived  measurements  of  intensity  and  phase  scintillation  (as  characterized  by  the  S4  and 
0+  indices  and  the  T  and  p  phase-spectrum  parameters)  it  is  necessary  to  convert  the 
measurements  into  estimates  of  the  height-integration  strength  of  the  electron-density 
irregularities  which  cause  the  scintillation,  denoted  C^L.  The  C*L  estimates  can  then  be  used  to 
modify  the  model  of  C*L  within  WBMOD  to  provide  an  improved  specification  of  ionospheric 
scintillation  for  AWS's  use.  The  source  of  data  will  be  the  modified  GPS  receivers  which  will 
make  up  the  AWS  TISS  planned  for  deployment  in  the  mid  1990s. 

Two  programs  have  been  developed  and  delivered  to  the  AWS  for  this  purpose.  Program 
GBCKL  (Ground-Based  C*L)  converts  TISS -derived  measurements  of  T,  p,  S4  and  ot  to  esti¬ 
mates  of  C*L  using  the  same  single-regime  power-law  propagation  theory  used  to  calculate  scin¬ 
tillation  effects  from  C*L  in  the  WBMOD  model  These  estimates  of  C*L,  with  estimates  gener¬ 
ated  from  in  situ  observations  of  the  irregularities  from  the  DMSP  Scintillation  Meter  instrument, 
are  used  in  program  IONSAM  to  produce  an  analysis  of  the  global  distribution  of  C*L  for  use  m 
the  WBMOD  scintillation  model  These  two  programs  are  part  of  the  Ionospheric  Scintillation 
Specification  and  Prediction  System,  which  has  been  tested  and  installed  at  the  USAF  Space 
Forecast  Center  (AFSFC)  at  Falcon  AFB,  CO.  Two  versions  of  the  IONSAM  code  were  de¬ 
livered,  one  that  works  within  the  DECWindows  windowing  environment  and  produces  several 
displays  and  plots  of  the  results  of  the  analysis,  and  one  with  a  simple  command-driven  interface 
which  is  capable  of  running  on  any  system  with  a  FORTRAN  77  compiler. 

The  software  was  formally  tested  at  the  AFSFC  and  the  results  of  the  testing  were  docu¬ 
mented  in  a  Test  Analysis  Results  document  which  was  provided  to  AWS.  The  run  times  of  the 
GBCKL  and  IONSAM  program  codes  extracted  from  the  test  runs  are  summarized  in  Table  2.  A 
TELSI  record  is  a  single  data  record  which  is  to  include  a  15-minute  summary  of  scintillation 
observations  from  a  single  TISS  site.  It  can  include  data  from  up  to  four  satellites.  An  SSIES 
record  is  a  C*L  estimate  generated  from  DMSP  Scintillation  Meter  observations.  During  the 
testing  period,  the  AFSFC  personnel  were  given  a  briefing  in  which  the  analysis  methods  imple¬ 
mented  in  GBCKL  and  IONSAM  were  described,  and  guidelines  for  using  the  software  to  meet 
their  needs  were  presented. 

Documentation  describing  the  use  and  structure  of  these  two  codes  was  also  delivered  to 
the  AFSFC.  Total  documentation  delivered  included  a  Test  Plan,  a  Test  Analysis  Report,  two 
User  Manuals  (one  for  each  program),  and  two  Program  Maintenance  Manuals  (one  for  each 
program).  Copies  of  the  documentation  were  sent  to  the  AFSFC  and  to  AWS/SCA  (Capt 
Jackson). 
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1.  Program  cbckl: 


17-station,  24-hour  run  (6,370  TELSI  records):  895s 
(2.2s  per  station-hour] 

2.  Program  ions  AM: 

Data  Input  Section: 

3-hour  TISS  (3,570  data  points)  and  SSIES  (1,153  data  points):  49s 
24-hour  TISS  (33,070  data  points)  and  SSIES  (8,642  data  points):  351s 

Analysis  Sections: 

Equatorial 

24-hour  data  set  (10,086  data  points):  2s 
Mid-latitude 

24-hour  data  set  (4,126  data  points):  Is 
High-latitude 

Startup:  2s  +  time  of  one  iteration  step 
Per  iteration  step: 

Non-binned  (4,904  data  points):  10s 
Binned  (323  data  bins):  0.8s 

Total  for  all  three  analysis  sections  for  a  typical  analysis:  63s 
24-hour  data  set,  equatorial  section  (10,000  data  points) 

24-hour  data  set,  mid-latitude  section  (4,000  data  points) 

3-hour  data  set,  high-latitude  section  (4,000  data  points) 
Non-binned  analysis,  5  iteration  steps 


Table  2.  AFSFC  testing  -  timing  results. 
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IV.  CONCLUSIONS 


A  large  set  of  TEC  and  scintillation  data  has  been  collected  from  mid-latitude,  sub-auroral, 
and  polar  stations  over  more  than  a  year  near  solar  maximum.  The  TEC  data  from  the  sub- 
auroral  station  (Lerwick,  Shetland  Island,  UK)  have  been  well  analyzed  and  are  shown  to  provide 
extremely  detailed  morphology  for  all  times,  seasons,  and  levels  of  magnetic  activity.  This  data 
set  includes  four  independent  measurements  in  four  different  directions  throughout  the  period, 
and  techniques  were  developed  to  display  these  data  to  show  the  perspective  of  one  monitor 
looking  nearly  overhead  with  two  other  monitors  located  north  and  south  of  the  station.  With  this 
technique,  latitudinal  variations  and  motion  of  the  ionospheric  trough  region  were  observed.  It 
was  seen  that  during  periods  of  activity,  in  most  instances  coincident  with  time  of  high  planetary 
Kp  index,  the  observations  are  consistent  with  equatarward  movement  of  the  irregularity 
boundary.  This  conclusion  supports  earlier  studies  performed  by  other  observers  (Tsunoda, 
1988).  Although  four-direction  data  were  not  available  at  the  other  sites,  other  data  sets  exist 
which  may  be  studied  in  conjunction  with  the  data  obtained  here.  For  example,  concurrent  130- 
and  400-Mhz  scintillation  data  were  obtained  at  Shetland  Island,  and  sounder  data  are  available  at 
Thule  and  Han  scorn  AFB. 
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APPENDIX  1 


ACTIVITIES  SUPPORTING  STATEMENT  OF  WORK  FOR  INVESTIGATION  OF 
HIGH  LATITUDE  IONOSPHERIC  STRUCTURES  NEAR  SOLAR  MAXIMUM 

During  the  22-month  study  period,  NWRA  maintained  and  operated  (or  directed  the 
operation  of)  data  collection  equipment  at  three  sites:  Thule  AB,  Greenland;  Lerwick  Obser¬ 
vatory,  Shetland  Island,  UK;  and  Han  scorn  AFB,  MA.  The  GFE  operated  and  serviced  by 
NWRA  consists  of  STEL-5010  single-channel  GPS  receivers,  a  TI-4100  four-channel  GPS 
receiver,  MX- 1502  Transit  receivers,  and  a  NIMS  code-free,  multi-channel,  GPS  receiver  manu¬ 
factured  by  AIR. 

The  STEL-5010  located  at  Shetland  Island  experienced  many  failures  during  this  study 
period.  Many  resources,  both  financial  and  manpower,  were  expended  to  maintain  the  site  in  an 
operating  condition.  This  impacted  the  quantity  and  quality  of  data  that  were  collected  at  the 
Shetland  Island  site  as  well  as  the  amount  of  data  from  all  sites  that  were  processed.  Emphasis 
was  placed  on  the  reduction  of  TI-4100  four  channel  GPS  Shetland  Island  data,  which  are  appli¬ 
cable  to  the  BMEWS  radar  at  Fylingdales,  UK. 

The  STEL-5010  single  channel  GPS  receiver  at  Han  scorn  failed  on  2  November  1991. 
Repair  was  attempted  using  parts  from  the  spare  STEL-5010  receiver.  The  spare  receiver,  having 
been  used  for  parts  several  times,  contained  no  usable  cards,  and  repair  of  the  Hanscom  receiver 
could  not  be  accomplished.  The  spare  receiver  was  sent  to  STEL  for  repair.  Several  months 
later,  a  repair  estimate  was  received  from  STEL  that  was  rejected  by  AFPL  as  too  costly.  In  view 
of  the  cost  vs.  benefits  analysis,  the  receiver  at  the  Hanscom  site  was  not  repaired. 

An  MX-1502  Transit  receiver  remained  in  operation  at  Hanscom  AFB.  This  receiver 
provides  latitudinal  relative  TEC  profiles  from  north-south  traveling  satellites  of  the  NNSS.  In 
order  to  determine  actual  TEC  levels,  these  data  require  that  complementary  absolute  TEC  data 
be  used  to  calibrate  the  relative  TEC  measurements,  as  was  done  for  UCW  Transit  data.  Data  are 
available  from  this  geographic  location  that  could  be  used  for  such  a  calibration  process.  From 
January  1991  through  June  1992,  data  from  a  geostationary  satellite  were  collected  using  a 
polari meter  located  at  Hamilton,  MA.  These  data  can  be  requested  through  Mr.  Jack  Klobuchar 
at  AFPL/GPIM,  Hanscom  AFB,  MA,  01731-3010.  Ionosonde  data  collected  at  Millstone  Hill, 
MA,  arc  available  in  several  different  formats  from  World  Data  Center,  Boulder,  CO,  and  could 
also  be  used  as  reference.  (Ionosonde  data,  collected  during  the  contract  period,  with  approxi¬ 
mately  90  percent  yearly  coverage,  is  also  available  from  Qanac,  Greenland,  and  Narssarssuaq, 
Greenland.) 

During  the  22-month  study  period,  NWRA  maintained  and  operated  (or  directed  the 
operation  of)  a  variety  of  satellite  receivers  at  Hanscom  AFB,  MA;  Lerwick  Observatory, 
Shetland  Island,  UK;  and  Thule  AB,  Greenland.  The  near  solar  maximum  database  in  Appendix  4 
details  the  availability  of  data  according  to  date,  site,  and  receiver  type,  and  indicates  the  level  to 
which  these  data  have  been  processed  and  the  quality  of  that  processed  data.  The  following  is  a 
chronology  of  the  problems  and  repairs  that  were  experienced. 
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22  December  1990  Shetland :  S TEL- 5010  receiver  failed.  Repaired  during  die  next  scheduled 
service  trip  on  19  February  1991. 

7*21  January  1991  Thule  :  Due  to  logistical  problems  associated  with  the  Gulf  War,  delivery 
of  data  tapes  to  Thule  ABwas  not  possible,  and  consequently  no  data  were  recorded  for  40  days 
during  the  period  11  January  1991  to  20  February  1991. 

18  January  1991  Hanscom:  Calibration  performed  on  the  STEL^5010  system. 

19  February  1991  Shetland  :  A  faulty  L2  microwave  oscillator  was  identified  as  the  source  of 
the  problem.  Replaced  with  a  nwrilar  oscillator  from  test  transmitter  until  arrival  of  a  permanent 
replacement.  Part  arrived  at  Shetland  on  5  March  91  and  was  installed.  New  5v  power  supply 
installed  in  GPS  receiver.  A  calibration  and  tune  up  was  successfully  performed  on  the  repaired 
receiver.  Regular  maintenance  performed  on  the  TI-4100  GPS  receiver. 

4  April  1991  Shetland  :  STEL-5010  receiver  failed. 

13>21  May  1991  Thule :  Service  trip  to  field  site.  Calibration  performed. 

26  May-2  June  1991  Shetland  :  Service  trip  to  field  site.  GPS  receiver  repaired  using  parts 
from  a  spare  receiver  sent  up  as  a  replacement  for  old  receiver.  Spare  receiver  inoperable  due  to 
rough  handling  during  shipping.  A  complete  calibration  of  the  STEL-5010  was  performed. 
Regular  maintenance  performed  on  the  Tl-4100  GPS  receiver. 

15  July  1991  Hanscom  :  Magnavox  MX- 1502  receiver  returned  to  AFPL  from  ARL-UT. 
Installed  for  test  period  of  two  months  before  being  deployed  to  Shetland  Island  during  next 
service  trip. 

20  July  1991  Shetland  :  A  loose  cable  connection  was  detected  when  data  from  site  were 
processed.  It  was  traced  and  temporarily  fixed  by  site  personnel  until  the  next  scheduled  main¬ 
tenance  trip.  Also,  the  system  printer  exhibited  intermittent  failure,  causing  the  entire  data  col¬ 
lection  system  to  halt.  Operator  intervention  was  necessary  to  restart  the  system  after  such  a 
failure. 

2  November  1991  Hanscom :  The  STEL-5010  receiver  failed.  The  problem  was  traced  to  the 
antenna  pre-amp  switch  that  was  found  to  be  half  full  of  water.  This  shorted  out  the  amplifiers  in 
the  pre-amp  and  caused  some  internal  problems  in  the  receiver.  Repair  was  deferred  until  the 
return  of  the  spare  STEL-5010  receiver  from  STEL. 

4-11  November  1991  Shetland  :  Service  Trip  to  field  site.  Equipment  serviced,  calibration 
performed,  and  the  MX- 1502  Transit  receiver  installed.  The  STEL  5010's  antenna  was  relocated 
from  the  roof  to  within  40  feet  of  the  TI-4100's  antenna  using  400  feet  of  low-loss  antenna  cable. 
Regular  maintenance  performed  on  the  TI-4100  GPS  receiver. 

8  February  1992  Shetland  :  The  TI-4100  receiver  was  turned  off  and  shipped  back  to  ARL- 
UT  for  calibration  and  reprogramming  before  deployment  to  Shemya,  AK,  on  another  effort 
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8-12  March  1992  Thule  :  Service  trip  to  field  site.  Calibration  performed.  Plexiglas  dome  on 
the  antenna  was  found  destroyed  and  a  temporary  fix  implemented  while  replacement  was  sought 

17  March  1992  Shetland  :  The  VT-103  video  display  on  the  GPS  receiving  system  failed 
Data  collection  continued  while  a  replacement  display  unit  was  sent 

02  April  1991  Hanscom  :  A  new  multichannel  GPS  receiver  from  AIR  was  received  for  test¬ 
ing.  It  was  found  to  have  several  defects  and  was  returned  to  the  manufacturer  for  repair. 

15  July  1992  Shetland  :  Data  collection  halted  due  to  a  dead  blower  motor  on  the  Cipher  tape 
drive.  Replacement  tape  drive  was  sent  but  system  was  still  inoperable  due  to  the  need  for  a 
complete  system  tune  up.  A  decision  was  made  at  AFPL/GPIA  to  shut  down  the  STEL-5010 
GPS  receiver  system  due  to  cost  and  time  involved  in  getting  the  system  b*ck  in  working  order. 
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APPENDIX  2 


QUALITY  ASSURANCE  EFFORTS  FOR  DATA  ACQUIRED  IN  THE 
INVESTIGATION  OF  HIGH  LATITUDE  IONOSPHERIC  STRUCTURES  NEAR 

SOLAR  MAXIMUM 

In  order  to  determine  TEC  values  accurately,  the  measurement  contributions  of  satellite 
biases  and  receiver  offsets  must  be  considered.  The  error  introduced  by  the  receiver  can,  in  the 
case  of  the  STEL-5010  single-channel  receiver,  be  determined  by  a  calibration  process.  This  is 
accomplished  by  injecting  a  known  signal  from  a  test  transmitter  directly  into  the  receiver  system 
at  the  pre-amp  and  measuring  the  response  in  volts.  Attenuators  are  inserted  between  the  test 
transmitter  and  the  pre-amp.  Attenuation  is  increased  until  the  receiver  begins  to  lose  lock  on  the 
weak  signal.  The  attenuation  level  is  backed  off  30  dB.  Several  steps  of  delay  are  then  intro¬ 
duced  and  the  corresponding  receiver  response  is  measured.  From  this  set  of  input  dB  of 
attenuation  and  output  voltages,  a  linear  relationship  is  determined.  This  calibration  is  used  in  the 
processing  of  the  raw  data  in  the  programs  GPSPCNW3.FOR  and  GPSPC3.FOR.  Each  recorded 
voltage  measurement  is  converted  to  dB  of  attenuation  according  to  the  linear  relationship 
determined  by  the  calibration.  The  process  in  repeated  for  both  the  LI  and  L2  channels. 

A  calibration  is  also  done  on  the  DGD  output  Signal  delay  is  introduced  in  6  NS  (INS  - 
2.852  TEC)  steps  through  a  range  of  278.83  TEC,  and  the  output  of  the  DGD  channel  is 
measured.  This  set  of  inputs  and  outputs  produces  a  best-fit  linear  relationship,  which  is  applied 
to  the  raw  data  measurements  by  the  programs  GPSPCNW3.FOR  and  GPSIK3.FOR  to  deter¬ 
mine  the  actual  values  of  DGD. 

In  addition  to  the  calibration  adjustments  to  the  raw  TEC  data,  an  adjustment  is  also  made 
to  correct  for  satellite  biases.  Satellite  biases  are  the  result  of  the  actual  delay  between  the  LI  and 
L2  signal  transmissions  from  the  satellite.  This  correction  is  applied  by  removing  the  bias  values 
from  slant  TEC  measurements  during  processing  through  GPSPCNW3.FOR  and  in  the  plotting 
programs  before  conversion  to  equivalent  vertical  TEC. 

Plotting  of  the  data  provides  a  quality  check  on  the  collection  equipment's  performance. 
Several  programs,  described  in  Appendix  3,  are  available  to  plot  the  data  as  24-hour  profiles,  24- 
hour  several-day  overplots,  and  individual  satellite  passes.  The  azimuth  and  elevation  of  the 
satellite  can  be  displayed  on  the  pass  file  plots,  while  satellite  latitude  and  longitude  can  be  shown 
with  the  24-hour  format  plots. 

While  processing  the  December  1991,  January  1992,  and  February  1992  data  from 
Shetland  Island,  it  was  noted  that  the  24-hour  TEC  profiles  exhibited  negative  values  and  dis¬ 
tortion  which  might  be  attributable  to  incorrect  values  of  satellite  biases  and  receiver  offset  being 
used  by  the  plotting  routines  Vertec3.for  (single-channel  data)  and  Vertec4.for  (4-channel  data). 

An  analysis  was  undertaken  to  find  values  for  satellite  biases  and  receiver  offset  that 
would  yield  TEC  curves  having  reasonable  shape  and  minimum  TEC  level.  First,  several  weeks 
of  single-channel  data  were  overplotted  in  one- week  sets.  The  assumption  was  made  that  receiver 
offset  and  satellite  biases  were  approximately  constant  over  this  period.  The  character  of  the  TEC 
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curves,  with  respect  to  expected  behavior,  was  examined.  Changes  to  receiver  offset  and 
individual  satellite  biases  were  estimated.  The  data  were  rerun  using  the  revised  values  and 
analyzed  again.  The  several  week  period  (three  months)  allowed  examination  of  most  satellites 
during  some  nighttime  hours,  due  to  the  satellites  appearing  approximately  four  minutes  earlier 
each  day.  The  conversion  to  equivalent  vertical  produced  distorted  data  when  the  combined 
satellite  bias  and  receiver  offset  were  incorrect,  guiding  the  revised  estimate  for  the  next  iteration. 
This  process  was  repeated  until  an  offset  and  a  set  of  bias  figures  for  the  block- 1  satellites  that 
resulted  in  reasonably  shaped  TEC  curves  were  determined. 

Four-channel  GPS  data  from  the  same  time  period  and  site  were  studied  next  Using  the 
block- 1  satellite  biases  determined  from  the  single-channel  data,  the  four-channel  data  were 
analyzed  to  select  a  receiver  offset  that  resulted  in  similar  TEC  profiles  for  the  same  satellites  at 
the  same  times  as  the  single  channel.  Once  this  was  done,  the  biases  for  block-2  satellites  were 
estimated  by  referencing  a  block-2  satellite  at  a  given  latitude,  longitude,  and  time  to  a  block- 1 
satellite  at  very  nearly  the  same  latitude,  longitude,  and  time.  The  data  were  reprocessed,  refining 
the  biases  each  time.  After  several  iterations,  the  estimated  values  of  receiver  offset  and  satellite 
biases  produced  24-hour  plots  of  TEC  whose  behavior  exhibited  expected  diurnal  and  seasonal 
variations,  no  negative  TEC  values,  and  credible  minimum  levels.  Table  3  lists  the  empirically 
derived  satellite  bias  values. 


Table  3. 

Derived  Satellite  Biases 


Satellite  No. 

Bias  Value  (TECu) 

Satellite  No. 

Bias  Value  (TECu) 

2 

2.20 

15 

1.96 

3 

10.67 

16 

7.77 

6 

7.33 

17 

2.75 

8 

0.00 

18 

-3.99 

9 

-0.97 

19 

7.76 

11 

18.01 

20 

2.59 

12 

10.41 

21 

17.54 

13 

10.73 

23 

3.60 

14 

2.04 

24 

5.00 

In  order  to  evaluate  the  accuracy  of  these  empirically  determined  offset  and  bias  values, 
several  single-channel  data  sets  were  overplotted  using  offset  values  which  varied  by  steps  of  +/-2 
TEC.  These  data  sets  are  shown  in  Figures  33  -35.  It  can  be  seen,  in  Figure  35  particularly,  that 
very  little  change  from  the  empirically  derived  bias  values  begins  to  introduce  distortion  in  the 
TEC  curves,  indicating  that  the  accuracy  of  these  values  is  about  ±  2  to  3  TEC  units. 

This  set  of  bias  values  determined  on  the  basis  of  December  91,  January  92,  and  February 
92  data  proved  to  be  valid  for  data  back  to  March  91  and  ahead  to  March  92.  It  should  be  noted 
that  the  four-channel  1991  data  showed  a  gradual,  overall  shift  in  level  of  the  TEC  curves.  This  is 
believed  to  be  caused  by  a  drift  in  the  TI-4100  receiver  due  to  aging  of  components.  The  receiver 
seemed  to  be  stable  over  certain  time  periods,  ranging  from  one  week  to  three  months,  ami  an 
appropriate  offset  was  determined  for  each  of  these  periods  using  the  procedures  described  above. 
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For  days  060  - 174, 1991,  the  offset  was  determined  to  be  -S.91  TECu;  days  17S  -  230, 1991,  the 
offset  was  -10.91;  days  231  -  237,  1991,  the  offset  was  -15.91  TECu;  and  for  days  247,  1991 
through  039,  1992,  the  offset  was  -20.91  TECu.  The  single-channel  receiver  offset  was 
determined  to  be  -16.0  TECu. 
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I  i  lie  (  U  .  T  .  ) 


TIME  ( U. T  .  ) 


°  0  3  6  9  12  15  18  Zl  24 

T  1  ME  (  U  .  r .  ) 


Figure  33.  Satellite  bias  accuracy  analysis  -  Day  01  1992. 
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APPENDIX  3 


DESCRIPTION  OF  SOFTWARE  FOR  REDUCTION  AND  ANALYSIS  OF  NEAR 
SOLAR  MAXIMUM  HIGH  LATITUDE  DATABASE 


The  software  developed  during  this  effort  reduces  the  raw  data  to  files  from  which  24- 
hour  plots,  several-day  overplots,  pass-file  plots,  latitude-separated  plots,  and  IPP-adjusted  time 
plots  of  equivalent  vertical  TEC  are  produced.  These  programs  can  be  used  with  any  GPS  data, 
such  as  that  collected  by  the  planned  TISS,  in  the  same  maimer  as  they  are  used  with  the  NWRA 
database:  to  test,  quality  check,  reduce,  analyze,  and  display  the  data. 

The  programs  developed  for  this  study  are  written  in  FORTRAN  for  use  with  Microsoft 
FORTRAN  for  PC  and  are  described  in  the  following  text  The  following  information  is 
organized  by  site  and  receiver  type.  Programs  to  process  data  collected  by  the  same  receiver  type 
at  different  sites  will  only  be  shown  with  an  underscored  title. 

A.  SHETLAND  ISLAND,  UK 

During  the  period  of  study,  three  receivers  were  operated  at  the  Shedand  Island  site:  an 
STEL-5010  single-channel  GPS  receiver,  a  TI-4100  four-channel  GPS  receiver,  and  an  MX-1502 
Transit  receiver. 

1.  STEL-5010  Single-Channel  GPS  Receiver 

The  STEL-5010  single-channel  GPS  receiver  collected  data  from  Block-1  GPS  satellites 
and  recorded  the  data  on  10-inch  magnetic  tapes.  Each  tape  holds  approximately  48  hours  of 
data.  DGD  between  LI  (1575.4  Mhz)  and  L2  (1227.6  Mhz),  DCP  between  those  two  signals, 
and  the  received  intensities  of  the  LI  and  L2  signals  were  recorded  at  20  Hz  as  binary  files.  The 
following  programs  are  used  to  process  these  data. 

a.  GPSPCNW3.FOR 

Input  to  this  program  is  the  raw  STEL-5010  GPS  data  from  the  10-inch  tapes.  This 
program  can  be  used  on  all  data  taken  with  the  STEL-5010  receiver  from  any  site  and  any  year. 
The  output  is  an  ASCII  file  containing  a  table  that  lists  universal  time  (seconds),  intensity  scin¬ 
tillation  (S4)  at  LI  and  L2 ,  slant  TEC,  rms  DCP,  azimuth  and  elevation  angle  of  the  satellite  from 
the  ground  station,  geomagnetic  latitude,  geomagnetic  longitude,  geomagnetic  time,  and  the 
penetration  angle  at  the  350  km  IPP.  A  header  at  the  beginning  of  the  file  lists  the  GPS  satellite 
ID  number,  Julian  date,  year,  and  the  name  of  the  receiving  site. 

The  data  are  one  minute  averaged  values  listed  every  20  seconds.  Geometry  data  and  the 
time  values  correspond  to  the  center  of  the  one-minute  averaging  interval.  The  RMS  DCP  data 
represent  a  measure  of  the  standard  deviation  in  dispersive  phase.  The  values  listed  are  RMS 
values  of  TEC  units  over  a  one  minute  averaging  interval  that  were  obtained  using  data  from  the 


55 


derampcd  Doppler  channel.  It  should  be  noted  that  the  deramped  Doppler  is  equivalent  to 
relative  TEC 

The  following  subroutines  are  used  and  are  described  as  follows: 

SUBROUTINE  RMSPHI 

Data  from  the  deramped  Doppler  channel  are  processed  to  obtain  one-minute  average 
RMS  values  which  are  written  to  an  intermediate  file. 

SUBROUTINE  NAV12 

Time  and  measurement  data  consisting  of  one-minute  averages  at  20  seconds  increments 
are  read  from  an  internal  file  and  stored  in  memory.  Keplerian  elements  applicable  to  die  GPS 
satellite  being  tracked,  measured  time  and  station  geodetic  coordinates  are  input  to  subroutine 
NAVSTAR.  NAVSTAR  returns  the  azimuth  and  elevation  of  the  satellite  at  the  ground  station 
and  the  geodetic  coordinates  of  the  satellite  latitude  and  longitude.  NAVSTAR  was  originally  a 
pass  planner  that  produced  files  containing  the  described  geometry  information  for  all  GPS 
satellites  observable  from  a  specified  ground  station  at  10  minute  increments.  The  modified 
NAVSTAR  subroutine  enables  the  user  to  directly  obtain  the  necessary  geometry  information 
without  having  to  access  external  files,  thereby  circumventing  an  additional  processing  procedure 
that  otherwise  would  be  critical. 

SUBROUTINE  MAGIPP(QNAME4) 

MAGIPP  is  a  geometry  conversion  program  that  inputs  receiver  coordinates,  measure¬ 
ment  time  (UT),  and  the  outputs  obtained  from  NAVSTAR  as  described  above,  and  produces  an 
intermediate  file  containing  penetration  angle,  geomagnetic  latitude  and  longitude,  and  geomag¬ 
netic  local  time  all  referenced  to  the  350  km  IPP.  In  accomplishing  this  task  an  additional  sub¬ 
routine,  IPP,  is  used  that  calculates  the  latitude  and  longitude  of  the  IPP,  given  die  receiver 
location  and  azimuth  and  elevation  look-angles  to  the  satellite. 

PROCESSING  ERROR  INFORMATION 

During  routine  data  processing,  errors  detected  are  written  to  disk  on  file  geomN_err 
where  N  is  incremented  for  each  run  (1-9  range) . 

Error  conditions  identified  are  as  follows: 

1)  almanac  data  necessary  for  updating  the  Keplarian  elements  are  identified 

2)  almanac  data  on  file  is  out-of-date  by  a  specified  number  of  days 

3)  almanac  data  file  does  not  match  the  requested  satellite  ID 

4)  elevation  angle  below  the  horizon 

5)  almanac  data  file  has  bad  format 

6)  failure  of  the  Kepler  subroutine  to  converge 
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b.  GPSPC3.FOR 

This  program  is  used  to  process  the  raw  STEL-5010GPS  data  stored  as  20  Hz  binary  data 
on  the  10-inch  tapes.  It  generates  four  ASCII  data  files  containing  DGD  and  DCP,  as  six-second 
averaged  values,  and  intensity  scintillation,  S4(L1)  and  S4(L2),  as  one-minute  averaged  values, 
calculated  every  20  seconds.  These  files  are  the  inputs  to  the  plotting  routines. 

c.  SHTREC2.FOR 

Shtrec2.for  phase  averages  the  DCP  data  using  the  DGD  as  reference.  Inputs  to  this 
program  are  the  DGD  and  DCP  files  created  by  GPSPC3.FOR.  Outputs  are  files  in  the  same  six- 
second  format  containing  phase  averaged  DCP  and  the  unaltered  DGD  data.  These  files  become 
the  inputs  to  the  plotting  programs. 

d.  VERTEC3.FOR 

Vertec3.for  produces  plots  of  the  DGD  and  DCP  data.  The  user  can  select  the  format  of 
the  plot  to  be  either  a  24-hour  equivalent  vertical  TEC  profile,  a  24-hour  overplot  of  several  days 
of  equivalent  vertical  TEC,  or  individual  pass  files  for  each  satellite  in  a  data  set  The  pass  file 
plots  can  display  slant  TEC  from  DGD,  slant  TEC  from  DCP,  and  vertical  TEC  from  phase-aver¬ 
aged  DCP  on  one  plot 

Inputs  to  Vertec3.for  are  the  two  files  generated  by  the  Shtrec2.for  program  and  a  file 
containing  satellite  ephemeris  data,  called  the  AZ-EL  file,  which  is  created  using  the 
RUNELEM.BAT  process.  The  output  is  a  24-hour  plot  or  pass-file  plots,  as  selected  by  die  user. 

For  Shedand  Island  S TEL-50 10  single-channel  data  only,  a  file  containing  a  table  of 
equivalent  vertical  TEC  and  IPP  latitude  and  longitude  at  five-minute  intervals  is  also  generated 
by  Vertec3.for.  These  tables  are  sent  to  UCW  for  use  in  calibrating  their  Transit  satellite  data.  In 
turn,  the  calibrated  Transit  measurements  become  the  input  to  the  Transit  Analysis  statistical 
evaluation  software  package,  described  below. 

e.  VERTEC4B.FOR 

Essentially  the  same  as  Vertec3.for,  this  program  can  also  plot  the  data  according  to  IPP 
adjusted  time,  as  opposed  to  UTC. 

f.  PLTS424.FOR 

This  program  plots  the  S4  files,  calculated  from  LI  and  L2,  generated  by  the 
GPSPC3.FOR  program  as  24-hour  plots  or  individual  satellite  pass  plots. 

g.  TRANSIT  ANALYSIS  SOFTWARE  PACKAGE 

RDP,  Inc.  of  Waltham,  MA,  has  created  software  to  evaluate  the  ETAC  ionospheric 

model  (based  on  the  Bent  ionospheric  model).  Transit  data  from  UCW  that  have  been  referenced 
to  the  GPS  data  collected  at  Shetland  Island  and  processed  by  NWRA,  are  compared  to  model 
predicted  values  and  a  statistical  analysis  is  performed.  This  software  was  modified  to  also  accept 
as  input  single  channel-format  GPS  absolute  TEC  data  from  the  TI-4100  receiver. 
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h.  RUNELEM.BAT 

The  purpose  of  this  batch  file  is  to  create  the  AZ-EL  file  containing  the  satellite  ephemeris 
data  needed  by  the  plotting  programs.  The  files  needed  by  this  process  are;  (1)  an  file 

from  the  Holloman  AFB  bulletin  board  covering  the  period  of  time  of  interest,  (2)  a  site  file 
containing  the  name,  latitude,  longitude,  elevation  of  the  site,  and  the  start/stop  Julian  dates  of  the 
desired  az-el  coverage,  and  (3)  the  programs  Gpselemlor  and  Navi  l.for.  The  output  is  an  az-el 
file  containing  the  ephemeris  data  needed  for  plotting. 

1.  GPSELEM.FOR 

Gpselem.for  reformats  the  almanac  downloaded  from  the  Holloman  AFB  bulletin  board  so 

that  it  can  be  used  by  the  Navi  l.for  program. 

2.  NAV11.FOR 

Navi  1. for  produces  an  az-el  file  that  contains  the  azimuth,  elevation,  latitude,  and  lon¬ 
gitude  of  each  satellite  at  10  minute  intervals,  for  each  day  in  the  period  specified  in  the 

site  file. 

L  PLTREC7C.FOR 

This  utility  program  allows  for  manual  correction  of  discontinuities  in  the  phase  data 
(DCP)  that  is  in  the  single-channel  (STEL  5010)  file  format  With  the  data  graphically  displayed 
on  the  screen,  the  operator  uses  the  cursor  to  mark  the  discontinuity.  The  program  will  remove 
the  discontinuity  by  connecting  the  discontinuous  data  and  changing  the  data  values  appropriately. 

2.  TI-4100  Four-Channel  GPS  Receiver 

A  TI-4100  four-channel  receiver  was  located  at  Shetland  Island  from  1  April  1990  to 
February  1992.  The  data  are  recorded  on  a  DC2000  40  MB  tape  cartridges  at  a  sampling  rate  of 
1  Hz.  Software  to  process  these  data  was  written  by  ARL-UT  for  the  CYBER  mainframe  com¬ 
puter,  which  was  decommissioned  in  August  1992.  Most  of  the  1991  data  were  converted  from 
original  format  to  a  single-channel  format,  i.e.,  the  STEL-5010  data  format,  so  that  they  could  be 
analyzed  and  plotted  in  the  same  manner  as  the  single-channel  data  for  ease  of  comparison. 
Because  the  DCP  data  experience  many  drop  outs  and  are  therefore  not  readily  suitable  for  the 
phase-averaging  process  without  discontinuity  correction  using  the  labor  intensive  Pitrec7c.for 
process,  the  DGD  data  only  are  plotted  as  equivalent  vertical  TEC. 

a.  SHTREC4.FOR 

This  program  is  similar  to  the  Shtrec2.for.  It  performs  the  phase  averaging  process  and 
generates  intermediate  files  which  are  the  inputs  to  Vertec4.for,  the  plotting  program.  If  the  DCP 
data  have  been  corrected  for  discontinuities  using  Pltrec7c.for,  then  equivalent  vertical  TEC  can 
be  generated  from  the  DCP.  If  discontinuity  correction  has  not  been  performed,  Vertec4.for  can 
produce  plots  of  equivalent  vertical  TEC  from  DGD  only.  Each  of  the  four  trackers  is  processed 
separately. 


b.  ERTEC4.FOR 

The  output  files  from  Shtrec4.for  become  the  inputs  to  Vertec4.for.  Also  needed  as  input 
is  the  AZ-EL  file  containing  the  satellite  ephemeris  data  from  the  time  period  of  the  data  to  be 
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plotted.  In  response  to  the  program's  questions,  the  user  determines  whether  the  output  is  a  plot 
of  equivalent  vertical  TEC  which  is  derived  from  DCP,  if  the  DCP  has  been  corrected  for  discon¬ 
tinuities  using  Pltrec7c.for,  or  a  plot  of  equivalent  vertical  TEC  from  DGD  only.  Twenty-four- 
hour  plots,  individual  pass  file  plots,  or  several  days  overplots  of  data  can  be  produced  in  order  to 
validate  the  data  or  study  TEC  behavior  over  a  specific  time  period.  Each  of  the  four  trackers  is 
processed  separately. 

A  file  containing  equivalent  vertical  TEC,  and  IPP  latitude  and  longitude  at  five  minute 
intervals  is  also  generated  by  Vertec4.for.  These  tables  are  sent  to  UCW  for  use  in  calibrating 
their  Transit  satellite  data. 


c.  SHETDIR4.FOR 

The  data  from  all  four  trackers  for  a  given  day  can  be  plotted  according  to  latitude  regions 
using  Shetdir4.for.  This  program  inputs  the  absolute  TEC  (DGD)  files  for  a  specified  day,  or 
series  of  days,  and  produces  24-hour  plots  of  the  data  from  the  selected  latitude  region;  north, 
overhead,  or  south  with  respect  to  the  receiving  station.  This  is  helpful  in  discriminating  between 
the  mid-latitude/trough  and  auroral/polar  cap  behavior.  Data  can  be  plotted  versus  universal 
coordinated  time  or  IPP  adjusted  time. 


3.  MX- 1502  Transit  Receiver 


a.  MX2IP.BAT 

ARL-UT  has  provided  a  software  package  to  process  the  data  recorded  by  the  MX- 1502 
receiver.  The  MX2IP.bat  file  reads  the  data  from  the  data  cassette  and  generates  an  ASCII  data 
file  in  ionospheric  parameter  format  containing  observing  station  information,  satellite  infor¬ 
mation,  time,  refraction  count,  change  in  TEC  from  last  time  to  current  time,  signal  strength, 
azimuth,  and  elevation.  The  data  can  be  plotted  using  the  Ipplot  program. 

B.  THULE  AB,  GREENLAND 

One  STEL-5010  single-channel  GPS  receiver  was  operated  at  the  Thule  AB,  Greenland. 

1.  STEL-5010  Single-Channel  GPS  Receiver 


a. 

b. 


GPSPCNW3.FOR 

GPSPC3.FOR 


c.  SHTREC2.FOR 

d.  VERTEC3.FOR 

e.  YERTEC4B.FQR 
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t  PL1S424.FQR 


g.  RWELEM.BAT 

1.  GPSELEM.FOR 

2.  HAVll-FQR 

h.  ELIREC7C.FQR 

C.  PHILLIPS  LABORATORY,  HANSCOM  AFB,  MA 

Three  receivers  were  in  operation  at  PL  during  the  period  of  study,  an  STEL-5010  single¬ 
channel  GPS  receiver,  an  MX- 1502  Transit  receiver,  and  an  NIMS  nulti-channel  GPS  receiver 
manufactured  by  AIR.  The  AIR  receiver  was  acquired  in  APRIL  1992  and  shortly  thereafter 
returned  to  the  manufacturer  for  repair.  In  a  letter  from  AIR  dated  27  July  1992,  it  was  reported 
that  the  receiver  fails  to  meet  specified  accuracy  levels  due  to  antenna  design  problems  which  they 
plan  to  correct  Therefore,  the  data  from  this  unit  is  necessarily  of  doubtfull  quality  due  to  this 
design  flaw. 


1.  STEL-5010Single  Channel  GPS  Receiver 


b.  GESPC3.FQR 

c.  SHEREC2.FQR 

d.  VERTEC3.FQR 

e.  YERTEC4BJQR 

f-  PLTS424.FOR 

g.  RUNELEMJIAT 

1  GPSELEM.FOR 

2.  NAV11.FOR 

h.  ELIREC7C.F0R 

2.  MX- 1502  Transit  Receiver 

a.  M2QIP.BAT 
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APPENDIX  4 


NEAR  SOLAR  MAXIMUM  HIGH  LATITUDE  DATABASE 

Data  were  collected  at  three  sites:  Thule  AB,  Greenland;  Lerwick  Observatory,  Shetland 
Island,  UK;  and  Hanscom  AFB,  MA.  The  latitude  and  longitude  of  each  site  is  given  in  the  table 
below. 

SITE  DESCRIPTION 

SHETLAND  ISLAND,  UK  Latitude  60.16  Longitude  -1.15 

THULE  AB,  GREENLAND  Latitude  76.55  Longitude  -68.67 

HANSCOM  AFB,  MA  Latitude  42.30  Longitude  -71.16 

The  STEL-5010  single  channel  GPS  data  from  all  three  sites  are  recorded  at  20  Hz  on  ten 
inch  magnetic  tape  in  binary  format  These  data  are  currently  being  backed  up  on  eight  mm  tapes. 
The  processed  data  are  ASCII  files  stored  on  computer  hard  disk.  Strip  charts  displaying  die  LI 
and  L2  intensities,  and  DGD  and  DCP  data  as  they  were  recorded  are  also  available. 

The  STEL-5010  raw  data  can  be  processed  in  two  ways.  The  data  can  be  processed  into 
an  ASCII  file  in  a  tabular  format  containing  one-minute  averaged  values  at  20  second  intervals  of 
intensity  scintillation  (S4)  for  LI  and  L2,  slant  TEC  from  DCP,  rms  DCP,  azimuth,  elevation,  and 
geomagnetic  latitude,  longitude,  time,  and  penetration  angle  at  the  IPP.  The  data  can  also  be 
reduced  into  four  individual  ASCII  files  containing  DGD  and  DCP,  both  as  six  second  averages, 
and  intensity  scintillation  (S4)  from  LI  and  L2,  these  being  one-minute  averaged  values  at  20- 
second  intervals.  Software  was  developed  to  plot  these  individual  files  for  analysis,  display,  and 
quality  checking  purposes. 

TI-4100  data  were  recorded  in  binary  format  on  40-megabyte  DC2000  mini-data 
cartridges.  The  processed  data  files  are  stored  as  ASCII  files  on  hard  disks. 

The  Transit  data  were  collected  on  Verbatim  DC30  data  cassettes  in  binary  format 
Processed  data  files  are  in  ASCII  on  hard  disk. 

The  data  collected  by  the  AIR  GPS  receiver  are  stored  in  the  receiver’s  RAM  and  archived 
once  a  week  to  files  on  hard  disk. 

Each  of  the  following  25  pages  shows  one  month  broken  down  into  days.  Specified  by 
site,  the  listing  includes  a  column  for  each  receiver  type.  An  X  indicates  that  raw  data  from  that 
site,  on  that  day,  collected  by  that  receiver,  are  available.  An  O  appears  if  no  data  exist.  In  a 
separate  column,  Y  for  yes  and  N  for  no  indicates  if  the  data  have  been  processed.  The  last 
column  shows  the  quality  of  the  data  that  have  been  processed;  0  if  data  not  processed,  1  for 
good,  2  for  questionable  possibly  recoverable,  and  3  for  bad. 

Note  that  UCW  Transit  data  from  1992  exist.  When  a  detailed  listing  of  that  information 
is  obtained,  it  will  be  included  in  the  updated  log 
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Shetland  I  Shetland  I  Shetland  iHanscom  Hanscom  Thule  Data  Processed  Quality  Of _ 

GPS  GPS  Transit  GPS'S  Transit  GPS  From  Sites _ Processed  Data 

STI  5010  TI-4100  UCW  PL  STI  NIMS  Magnavox  STI  5010  1 1 2 1 3 1 4 1 5 1 6|  7 1 8 1  1  =  OK/2  =  7/3  =»  BAD 
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